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Abstract 

Current theory and experimental data on propellant 
utilization in electron-bombardment ion thrusters are 
reviewed. Because the majority of investigations have 
been conducted with mercury, the presentation empha- 
sizes that propellant. The results are presented in as 
general a form as possible to facilitate use in areas 
other than space propulsion. 

Introduction 


Electron-bombardment ion sources have been used 
in ion thrusters since 1960. 1 The need in electric 
space propulsion is for large ion currents, but with low 
enough current densities to permit the integration of the 
source with long-lived accelerators. The typical range 
of interest for exhaust velocities in ion thrusters is 
20 000 to 40 000 m/sec, which is low compared with 
most physics-experiment ion-beam applications. This 
exhaust -velocity range of interest results from a mis- 
sion compromise between high propellant mass at low 
exhaust velocities and high powerplant mass (usually 
solar-cell arrays) for high exhaust velocity. 

Neutral propellant atoms also escape through the 
accelerator during normal operation, and are the loss 
of particular interest in this investigation. Thruster 
efficiency is defined 

Vt = 2mP 


where (in SI or mks units) F is the thrust in newtons, 
m is the total propellant flow rate in kg/sec, and P is 
the total thruster input power in watts. In terms of the 
propellant flow that is ionized and ejected from the 
thruster, m^, and the propellant utilization 77 u (77 u = 
mj/m), Eq. (1) can be rewritten 


quF 2 

Vt 2rh 1 P 


( 2 ) 


The power efficiency is defined 

F 2 

2rh i P ’ 


(3) 


so that the thruster efficiency can also be written 




(4) 


The thruster efficiency is thus proportional to propel- 
lant utilization, and an improvement in propellant 
utilization can be translated directly into improved 
thruster performance. 

The purpose of this paper is to describe the cur- 
rent state of the art for propellant utilization in 
bombardment-thruster ion sources, to present experi- 
mental performance in as generalized and useful a 
form as possible, and to indicate where future improve- 
ments might be expected. Because the majority of 
experimental and theoretical investigations have been 
conducted with mercury, the performance will of 
necessity be slanted toward that propellant. The use 
of other propellants will also be discussed to assist 
the reader whose primary interest is other than space 
propulsion. 


Typical Operation 

During typical operation with mercury the neutral 
density corresponds to about 1(T 4 torr, or ~10 18 mer- 
cury atoms per m 3 at ~500° K ion-chamber wall tem- 
perature. The ions are only about 10 percent as dense 
as the neutrals (^10 17 /m 3 ) and the electron density 
equals that of the ions. All of these densities are 
approximate mean values for the ion chamber as a 
whole when operated at design conditions. Special 
locations, such as the vicinity of a hollow cathode could 
result in large departures from the values given. The 
ion-chamber electron distribution is usually treated 
theoretically as a monoenergetic "primary 1 1 popula- 
tion superimposed on a low energy Maxwellian back- 
ground population. The primary electrons actually 
have, of course, a spread in energy. The use of the 
monoenergetic approximation, though, greatly reduces 
the complexity of calculations for a two -group distri- 
bution while introducing no substantial errors. Roughly 
10 percent of the electrons are energetic primaries at 
normal operating conditions. The mean energy of, 
these primaries corresponds to about the total potential 
difference across the discharge. The background 
Maxwellian population has a temperature of about 5 eV. 

The physical cause for a two -group electron dis- 
tribution resides in the collision probabilities for pri- 
mary electrons emitted by the cathode. For the typical 
30 to 40 eV primary electron, the probability of an 
ionizing or exciting collision is about 10 times as great 
as the probability of Coulomb randomization with the 
low energy background electrons. The high-energy 
primary population therefore comes from cathode 


TM X 07921 



emission and is decoupled from the low-energy 
Maxwellian population. 

Roughly half of the ions produced within the ion 
chamber come from primaries, while the rest come 
from the ’’tail" of the low-energy Maxwellians. 2 At 
high utilization conditions, the fraction of primaries 
can increase substantially, resulting in most of the ions 
coming from primaries at these conditions. 3 

The Debye shielding distance is ~5xl0~ 3 cm for' 
typical ion-chamber conditions. This small Debye 
shielding distance, when compared with ion-chamber 
dimensions, justifies the assumption of equal electron 
and ion densities. The bulk of the ion-chamber plasma 
is within a few volts of a uniform potential. The level 
of this nearly uniform potential is typically ^5 V above 
anode potential. Because of the small Debye shielding 
distance, sheaths near the anode and cathode are quite 
thin. 

The ion density is normally much less than the 
neutral density, but far more ions than neutrals leave 
through the accelerator system. This apparent dis- 
crepancy is because of the difference in velocities 
within the ion chamber for these two species. The 
neutral velocity is controlled by the wall temperature 
of ~500° K. The ion velocity results from electro- 
static fields within the ion-chamber plasma and is the 
minimum value found necessary by Bohm 4 for sheath 
stability at boundaries of ion production regions. 

The rapid escape of ions is also involved in the 
production of multiply charged ions. Multiple ioniza- 
tions, if significant, are produced by the initial col- 
lisions of electrons with neutrals. That is, ionizing 
collisions of electrons with previously ionized atoms 
are infrequent and can generally be ignored. Because 
a variety of charge -to -mass ratios produces a spread 
in ion beam velocity and a loss in power efficiency, it 
is customary to operate ion-thruster discharges at 
voltages such that only a small fraction of multiply 
charged ions are produced. 

The magnetic field strengths employed vary in- 
versely with thruster diameter. 5 For a 10-cm diam 
anode, the typical field strength is 1 to 5xl0” 3 T 
(10 - 50 G). This field strength is sufficient to contain 
the energetic primary electrons, but has little effect 
on the ion trajectories. 

Electron-bombardment thrusters have been oper- 
ated on propellants other than mercury. Next to 
mercury, cesium has been the most used propellant. 6 * 7 
A much smaller amount of operation has been obtained 
with gaseous propellants. An early study was made 
with gaseous propellants in a 10-cm diam thruster. 8 
More recently a study has been made using a divergent- 
field thruster. 9 The operation with cesium departs 
sharply from the preceding description for mercury. 


The lower discharge potential differences with cesium 
(~10 V) result in rapid randomization of electrons 
emitted by the cathode. As a result, a single 
Maxwellian distribution is found for the electrons 
instead of the two -group distribution obtained with mer- 
cury and most other propellants. The Coulomb ran- 
domization of electrons varies inversely as the square 
of electron energy. The reduced discharge voltage of 
cesium relative to mercury is therefore the cause of 
this electron-distribution difference. For discharges 
over about 20 V (of interest for most propellants) the 
assumption of the two -group electron distribution should 
be reasonable. 

The temperature of the Maxwellian background 
electrons is related to the excitation cross section. In 
mercury the excitation cross section becomes signifi- 
cant near 5 eV. It is not surprising, then, that meas- 
ured Maxwellian temperatures are typically 5±2 eV. 
Estimates for Maxwellian temperatures with other pro- 
pellants should therefore be guided by excitation cross- 
section knowledge — if available. 

Primary-Electron Region 

Recent high-performance ion-chamber designs have 
used sharply divergent field shapes in which the field 
strength passes through zero at an axial location near 
the accelerator and screen grids. Examination of ion- 
chamber probe surveys 2 * 10 have shown that the primary 
electrons are not found throughout the ion chamber. 
Instead, the primary electrons are concentrated in a 
region defined approximately by magnetic -field lines 
(guiding centers) that do not intercept the anode. In a 
divergent-field thruster this primary-electron region 
is located close to the accelerator and screen grids, as 
indicated by the cross-hatched region in Fig. 1. In this 
region are the two subregions, one that can be reached 
directly by electrons emitted by the cathode and one that 
can be reached only after collisions. The Maxwellian 
electrons have much less energy than the primary elec- 
trons, and thus can be largely reflected by the anode 
sheath. The Maxwellian electrons are therefore not 
confined to field lines that do not intercept the anode, 
and occupy the entire ion-chamber volume. 

The overriding importance of the primary-electron 
region to ion-chamber performance has been shown by 
the wide range of configuration changes that can be made 
outside this region with only small effects on perform- 
ance. 3 

An early attempt to predict bombardment-thruster 
utilization assumed uniform conditions throughout the 
ion chamber. 14 Such an approach is clearly not suitable 
for divergent -field designs. A more recent study of 
propellant utilization was concerned with conditions in 
only the primary -electron region. 3 The theory pre- 
sented in this recent study is summarized in the next 
section. 
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Figure L - Primary-electron region (shown by crosshatching) 
of divergent-field thruster. 


This equation can be used in a simplified form by as- 
suming a limiting condition where all electrons are 
primary electrons. Plasma neutrality then requires 
that ni /n p = 1 . To evaluate vj /v p we use the mini- 
mum ion velocity from the Bohm criterio n. This mini- 
mum ion velocity can be written Vi = v p /m e /2mj, 
where v p is the mean primary -electron velocity, and 
m e and m$ are the electron and ion masses. With 
these substitutions in Eq. (7), 


V m e /2m t 
^(Vp /Ap) 


( 8 ) 


All parameters of Eq. (8) are constants, giving a con- 
st ant value for neutral density iIq at maximum utiliza- 
tion. (Note that the ratio V p /A p is a characteristic 
length for the primary-electron region, and that this 
length should vary in a roughly proportional manner 
with ion-chamber diameter.) A constant neutral den- 
sity upstream of the accelerator system implies a con- 
stant loss of neutrals at maximum propellant utilization 
regardless of the total mass flow. 


A more sophisticated relationship can be obtained 
by using the previously mentioned combination of mono- 
energetic primaries and background Maxwellians for 
the electron distribution. A modified Bohm criterion 
has been derived by Masek 2 based on this two -group 
electron distribution. Using this approach, Eq. (7) 
becomes 


Maximum-Utilization Theory 


The density of primary electrons increases with 
discharge current. 3 For a maximum-utilization condi- 
tion at high discharge current, then, ion production 
should be approximated by considering only primary 
electrons. Using mean values for the primary-electron 
region, the total ion production rate N from primary 
elections is 


N = n pVpn 0 o-Vp, (5) 

where np and v p are the primary electron density and 
velocity, and a are the neutral density and ioniza- 
tion cross section, and V p is the volume of the pri- 
mary electron region. For the total loss rate of ions 
from this region, 


_ V m e /2m"i T m (1 +n p /n m ) 3/2 
°-a(V p /A p , JT p n p /n m ’ U 

where T m and are the Maxwellian electron tem- 
perature and primary electron energies — both in eV, 
and np/n m is the primary-to -Maxwellian density 
ratio . 

To compare the results of Eqs. (8) and (9), typical 
values can be used for the various constants. The mass 
of a mercury ion can be used for m^, 5X10 m^ used 

for a, and 40 eV and 5 eV used for and T_ 

_ * , . , . P m* 

These values yield 


2.34x1016 

° = V*P 


( 10 ) 


N = n^Ap, (6) 

where and are the ion density and velocity 
toward the outer boundary of the primary -electron 
regions and A p is the area of that outer boundary. 
Equating of the loss and production rates followed by 
solution for neutral density yields 

(ni/npXvj/vp) 

n ° = “olVp/ApT- (7) 


from Eq. (8) and 

0.827X10 16 (1 +n p/ n m) 3/2 
"° = V p /Ap n p /„ m (U) 

from Eq. (9). Equations (10) and (11) can be multiplied 
by Vp/A p . The i*esultof this operation for Eq. (11) is 
plotted in Fig. 2, while the result for Eq. (1.0) is indi- 
cated on the ordinate of Fig. 2. The value of 
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PRIMARY -TO-MAXWELLIAN DENSITY RATIO, n p /n m 


Figure 2. - Variation of neutral -density parameter with 
primary-to -Maxwellian density ratio from equation (11). 
Primary energy, 40 eV; Maxwellian temperature, 5 eV. 

no(V p /Ap) from Eq. (11) is seen to be higher than the 
same parameter from Eq. (10) for all values of 
n p/n m < I- Equation (11) is felt to be a more accurate 
representation of conditions in the primary electron 
region because of a more accurate electron-distribution 
representation. On the other hand, Eq. (10) requires 
less knowledge of ion-chamber conditions, and is there- 
fore easier to use. 

Propellant-Utilization Data 

Typical ion chamber performance is shown in 
Fig. 3. The curve for each propellant flow rate is rela- 
tively flat below a "knee." Above this knee the curve 
rapidly becomes nearly vertical. The location of the 
knee and the nearly vertical portion above the knee de- 
pends on propellant flow rate. 



Figure 3. - Typical ion-chamber performance data fora 
bombardment thruster (from ref. 3). 


The investigation that resulted in Eqs. (8) and (9) 
also compared these equations with experimental ion- 
chamber data. 3 The ion-chamber that was used for 
these experimental performance data is shown in Fig. 4. 
(The primary-electron region for this ion chamber was 
shown in Fig. 1. ) Important features of the design were 
movable and removable parts that permitted different 
ion-chamber lengths and diameters without affecting the 
primary-electron region. Also, a refractory metal 
cathode (tungsten wire) was used to avoid either the 
evaporation of oxide that can occur at high emission or 
the propellant concentration that results from the flow 
through a hollow cathode. 

Maximum-utilization data obtained from the ion- 
chamber of Fig. 4 is shown in Fig. 5. The definition 
of maximum utilization used for these data was the 



Figure 4. - Ten -centimeter thruster with adjustable ion chamber. 


utilization at 1000 eV/ion. This level of discharge loss 
put the performance well into the near-vertical portion 
of the ion-chamber performance curves — as indicated 
by Fig. 3. The level of neutral propellant loss obtained 
(Fig. 5(b)) was nearly a constant for a wide range of 
ion-chamber configurations and propellant flow rates. 
The level of the neutral loss could be made to agree 
with both Eqs. (8) and (9) by proper selection of con- 
stants. Equation (9) is more complete, and should be 
more suited to extrapolation to new operating conditions, 
but it requires knowledge of T m , ,£p, and n p /n m that 
is usually available only from Langmuir probes. 

The data of Fig. 5 are quite interesting and 
instructive in that they show the constant neutral loss 
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(B) NEUTRAL LOSS. 


Figure 5. - Correlation of utilization data from 
10-cm ion chamber. Discharge loss is 
1000 eV/ ion. Different symbols indicate dif- 
ferent shapes outside primary-electron 
region. 


predicted by Eqs. (8) and (9). Practical applications, 
though, center around the "knee” of the discharge -loss 
curve (Fig. 3). If one operates an ion chamber at a 
utilization below that of the knee, then the propellant 
flow rate could be reduced with only a small increase 
in discharge loss. Conversely, if one is operating at 
some high-discharge loss above the knee, then the dis- 
charge loss could be reduced with only a small decrease 
in utilization. 

Accordingly, available ion thruster data were ex- 
amined with a view toward finding data obtained: 

(1) near the discharge-loss knee, (2) under good enough 
experimental conditions (particularly facility size) to 
assure accurate propellant utilization measurements, 

(3) over a wide range of propellant flow rate, and 

(4) with strongly divergent magnetic -field designs. In 
addition to the 10-cm ion-chamber design shown in 
Figs. 1 and 4, data were used from chambers with 
5-cm, 15-cm, and 30-cm diameters. All three of 
these ion-chamber designs were the result of much ex- 
perimental optimization, and all used hollow cathodes. 
The primary-electron regions an,d general ion-chamber 
shapes for these three thrusters are shown in Fig. 6. 

All of these ion chambers (including the one shown in 
Fig. 4) were used with about 50-percent open accel- 
erator grids and 65- to 70 -percent open screen grids. 
The effective sharp -edged orifice area for all of these 
accelerator systems was estimated at one-quarter of 
the beam area. More detailed descriptions of the 

5-cm, 15-cm, and 30-cm thrusters can be found in 

12 

another paper presented at this conference. 
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(A) ANODE DIAMETER, 5 CM; V p /Ap, 0. 7 CM. 




(C) ANODE DIAMETER, 30 CM; V p /Ap, 2. 5 CM 

Figure 6. - Primary-electron regions (shown by crosshatching) of current divergent- 
field bombardment thrusters. 


Data obtained near the knee of the performance 
curves are shown in Fig. 7 for the four ion cham- 
bers. 13 ^ 6 With the exception of the 5-cm data, all 
data were obtained with discharge losses of 250-300 
eV/ion. The 5-cm data were obtained at 570 eV/ion. 

The 5-cm design represents a compromise toward 
improving the low utilization associated with small size 
at the expense of a higher discharge loss. This com- 
promise is evident in the deeper (in terms of length-to- 
diameter ratio) primary-electron region of the 5-cm 
design shown in Fig. 6. With this deeper region the 
value of Vp /A p is increased, but at the expense of a 
larger fraction of ions produced going to the ion- 
chamber walls instead of being accelerated into the beam. 
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(A) PROPELLANT UTILIZATION. 


where K includes all constants that are not a function of 
the ion source configuration or the propellant. With the 
substitution of Eq. (8) for r^, Eq. (13) becomes 


KA o y m e /2m i 
m[ /2 <7(Vp/A p ) 


(Note that for all practical purposes ) Finally, 

by redefinition of K, 



.01 .1 1 10 
PROPELLANT FLOW RATE, A -EQUIVALENT 

(B) NEUTRAL LOSS. 


N 0 = KAo/ mi cr(V p /Ap). (15) 

When solved for K, 

K = N o m i 0 ’(V p /A p )/A o . (16) 

Thus the product Noin^^p /A p )/A 0 should be nearly 
constant for all ion chambers. 

The data of Fig. 7 were replotted in Fig. 8 in terms 
of the neutral loss parameter N 0 m^cr(V p /Ap)/Ao. For 
convenience, the neutral loss rate N 0 was expressed 
as equivalent amperes (the current that would be meas- 
ured if each atom had one electronic charge) , and m* 
was expressed in atomic mass units. The remaining 
parameters cr, (Vp/A^), and A c were expressed in 
m 2 , m, and m 2 . 


Figure 7. - Utilization data for the "knee” of the discharge -loss curve. 

Mercury propellant. 

The utilization efficiency (Fig. 7(a)) clearly de- 
creases with total propellant flow rate for each of the 
ion chambers. The corresponding loss in neutrals is 
shown in Fig. 7(b). The neutral loss wanders over a 
small range for each thruster, sometimes increasing 
and sometimes decreasing with increasing total propel- 
lant flow rate. In going from the 1000 eV/ion data of 
Fig. 5 to the "knee" data of Fig. 7, the increase in 
practical interest has clearly been accompanied by an 
increase in scatter for each thruster. The data of 
Fig. 7, though, still indicate that neutral loss is ap- 
proximately constant for the "knee" data. 

A suitable parameter was desired to facilitate 
comparison of operation with different propellants. 
Equation (9) would have been used, but it requires 
more information than is generally available. Instead, 
Eq. (8) was selected to be the basis for deriving this 
parameter. The neutral density of Eq. (8) can be ex- 
pressed as a neutral propellant loss rate N Q when 
combined with an effective sharp-edged orifice area A 0 
for the accelerator system a nd the aver age magnitude 
of neutral velocity (v^ = Y 8 kT/irm). 

K = V 0 A 0 /4 (12) 

The wall temperature should be roughly the same for 
all ion chambers. We should therefore be able to write 

N 0 = Kn 0 A 0 /m 1 / 2 (13) 


In addition to the mercury propellant data of Fig. 7 , 
other gas data are plotted for the 15-cm ion chamber.^ 
Xenon and krypton were used in an unmodified 15-cm 
thruster, with knee data again obtained in the 250 - 
300 eV/ion range. A modified 15-cm ion chamber was 
used with argon. This modification was to block screen 
holes except for the center 7 . 5-cm diameter of the 
accelerator system. Inasmuch as this modification 
also blocked much of the ion extraction, the discharge 
loss was increased to about 750 eV/ion. Values used 
for the ionization cross section in the neutral-loss 
parameter were 5x10 “ 20 m 2 for mercury, 5x1 0~ 20 m 2 
for xenon, 4. 5xl0” 20 m 2 for krypton, and 3xl0“ 20 m 2 
for argon. These values were selected to be near the 
maxima, but still at low enough electron energies that 
double ionization would be negligible. 17 These gas 
propellant data were included in Fig. 8. 
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PROPELLANT FLOW RATE, A EQUIVALENT 


Figure 8. - Neutral -loss parameter for the "knee" of the discharge -loss 
curve. 
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Because the mercury data used corresponded to 
considerable thruster optimization, other gas propel- 
lant data were selected for good propellant distribution, 
near-optimum magnetic -field strength, and a moderate 
discharge voltage (to minimize multiple ionization). 
Further, no data were used for propellants that re- 
quired much higher discharge losses (neon and helium), 
because it was felt that the higher losses were indica- 
tive of some other than normal mode of operation. (It 
appeared that the neutral flow required for the hollow 
cathode to operate with neon and helium was too much 
for good ion-chamber performance. This viewpoint is 
supported by the fact that the earlier gas-thruster in- 
vestigation Actually gave better performance with these 
gases. 8 ) 

One can conclude from Fig. 8 that for the range of 
propellant flow rates investigated, the range of propel- 
lants used, and a strongly divergent magnetic -field 
design, the neutral loss parameter is about 
2. 1±1. Ixl0“ 18 . The degree of correlation appears 
good enough to permit use of this parameter as an ap- 
proximate design tool. Because the data of Fig. 8 were 
obtained with at least some selection of operating con- 
ditions and/or optimization, the use of the neutral-loss 
parameter should be considered as merely an adjunct to 
good design. 

Earlier bombardment ion-source data using gas- 
eous propellants and only a moderate degree of 
magnetic -field divergence 8 were not included in Fig. 8. 
The magnetic field at the middle of the accelerator was 
about 80 percent of the maximum value in the chamber 
for this earlier data, while the field strength drops to 
nearly zero at the same location in strongly divergent 
designs. Some checks were made of this earlier data. 
The neutral-loss parameters were found to be two to 
three times as large as for strongly divergent designs. 

Concluding Remarks 

The primary-electron region plays an important 
role in the ionization processes of a bombardment- 
thruster ion chamber. From both analyses of the 
primary -electron region and experimental data, neutral 
propellant losses have been shown to be nearly constant 
at high discharge-loss conditions — regardless of total 
propellant flow rate. 

A region of more practical interest is the knee of 
the discharge loss curve, which generally represents 
the best compromise between discharge losses and 
neutral propellant losses. The neutral losses at the 
knee of the performance curve are not as constant for 
a given thruster as at the high discharge loss condition. 
But the neutral losses are still near enough constant to 
use for an estimate of utilization over a range of total 
propellant flow rate. 

A neutral-loss parameter was obtained. This 
parameter was experimentally evaluated for several 


different propellants and ion-chamber sizes that have 
been used with strongly divergent magnetic -field de- 
signs. The correlation indicates that this neutral-loss 
parameter may be useful as a design tool. 

Future improvements in performance of bombard- 
ment ion thrusters could come in part from proper 
shaping of the primary-electron region. A roughly 
hemispherical shape would appear near optimum; that 
is, optimum from the viewpoint of maximizing the 
ratio Vp/Ap without excessively increasing Ap, as 
well as providing a fairly flat ion-beam profile. 

The shape of the primary-electron region shown in 
Fig. 1 would therefore appear desirable. Unfortunately, 
the field lines within this region result in a peaked ion- 
beam profile when the magnetic field is increased to 
minimize discharge losses. What is really needed, 
then, is a roughly hemispherical shape with near zero 
magnetic -field strength inside the primary-electron 
region and a strong magnetic field outside. Although 
these requirements may sound mutually exclusive, 
there is a magnetic -field design that may come close 
to satisfying them. This design is the multipole 
"picket fence" concept of plasma containment from 
fusion research. The multipole design has been used 
in a cesium thruster with high performance, 9 but the 
results with mercury have been mixed. 18 » 19 Some 
development of the multipole design appears likely for 
future high-performance bombardment thrusters. 
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